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SCATTERING OF COSMIC RAYS BY MAGNETOHYDRODYNAMIC 

INTERSTELLAR TURBULENCE 
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Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for substantial 
revisions in our understanding of cosmic ray transport. In this paper we use gyroresonance recently 
obtained scaling laws for MHD modes to calculate the scattering frequency for cosmic rays in 
the ISM. We consider gyroresonance with MHD modes (Alfvenic, slow and fast) and transit-time 
damping (TTD) by fast modes. We conclude that the gyroresonance with fast modes is the dominant 
contribution to cosmic ray scattering for the typical interstellar conditions. In contrast to earlier 
studies, we find that Alfvenic and slow modes are inefficient because they are far from isotropy 
usually assumed. 

PACS numbers: 98.70.Sa,52.30.-q, 52.35.Ra, 95.30.Qd 



I. INTRODUCTION 

The propagation of cosmic rays (CRs) is affected by 
their interaction with magnetic field. This field is turbu- 
lent and therefore, the resonant interaction of cosmic rays 
with MHD turbulence has been discussed by many au- 
thors as the principal mechanism to scatter and isotropize 
cosmic rays (§). Although cosmic ray diffusion can hap- 
pen while cosmic rays follow wandering magnetic fields 
H , the acceleration of cosmic rays requires efficient scat- 
tering. 

While most investigations are restricted to Alfven 
waves propagating along an external magnetic field (the 
so-called slab model of Alfvenic turbulence), obliquely 
propagating MHD waves have been included in H and 
later studies |l| The problem, however, is that 

the Alfvenic turbulence considered in their studies is 
isotropic turbulence, and this is contrary to the mod- 
ern understanding of MHD turbulence (||, see for a 
review and references therein). 

A recent study || found a strong dependence of scat- 
tering on turbulence anisotropy. Therefore the calcu- 
lations of CR scattering must be done using a real- 
istic MHD turbulence model. An important attempt 
in this direction was carried out in B. There Alfven 
modes were treated in the spirit of Goldreich-Shridhar 
(1995, henceforth GS95) model of incompressible tur- 
bulence and marginal scattering was obtained. However, 
a more accurate description is now available (|l(J) and 
thus there is a need to revisit the problem. Moreover, || 
did not consider compressible modes, while we show be- 
low that these modes provide the dominant contribution 
to the scattering. 



II. MHD STATISTICS 

MHD perturbations can be decomposed into Alfvenic, 
slow and fast modes (see PJ). Alfvenic turbulence is 
considered by many authors as the default model of in- 



terstellar magnetic turbulence. This is partially moti- 
vated by the fact that unlike compressible modes, the 
Alfven ones are essentially free of damping in fully ion- 
ized medium (see Jll], |l2|). 

Unlike hydrodynamic turbulence, Alfvenic turbulence 
is anisotropic, with eddies elongated along the magnetic 
field. This happens because it is easier to mix the mag- 
netic field lines perpendicular to the direction of the mag- 
netic field rather than to bend them. The GS95 model 
describes incompressible Alfvenic turbulence, which for- 



mally means that plasma (3 = P ga s/Pn 



infinity. It was first conjectured in ]13[ that GS95 scal- 
ing should be approximately true for moderately com- 
pressible plasma. For low plasma Cho & Lazarian |l4| ] 
(henceforth CL02) showed that the coupling of Alfvenic 
and compressible modes is weak and that the Alfvenic 
modes follow the GS95 spectrum |l5) . This is consistent 
with the analysis of observational data (EJ, 17 ). In what 
follows, we consider both Alfven modes and compressible 
modes and use the description of those modes obtained 
in CL02 to study CR scattering by MHD turbulence in 
a medium with energy injection scale L = lOOpc, density 
n = lfU 4 cm~ 3 , temperature T = 2 x 10 6 K. Recent ob- 



servations 1 18 suggest that matter in the galactic halos 
is magnetic-dominant, corresponding to low /? medium, 
here we choose (3 ~ 0.1. The injection length scale is im- 
portant as Alfvenic turbulence exhibits scale-dependent 
anisotropy that increases with the decrease of the scale. 

We describe MHD turbulence statistics by correlation 
functions. Using the notations from ||, we get the ex- 
pressions for the correlation tensors in Fourier space 



< B^B^k') > IBI = 5(k - k')M i3 (k), 
< Vi(k)B*(k') > /V A B = S(k - k')C (k), 
< Vi(k)v*{k') > fV% = S(k - k')ify(k), 



(1) 



where -B Q ,/3 is the magnetic field fluctuations. 

The isotropic tensor usually used in the literature is 



K tj (k) = C {S tJ - kikj/k 2 }k- n / 3 , 



(2) 



2 



The normalization constant Cq can be obtained if 
the energy input at the scale L is defined. As- 
suming equipartition, the kinetic energy density e/. = 
/ * 3 Eti K iiP Vl/2 ~ B 2 /8n, we get C = L- 2 / 3 /12tt. 
The analytical fit to the anisotropic tensor for Alfven 
modes, obtained in Jic| is, 

tf«(k) = ^^^^/^^(-L^^ii/fc^ 3 ), (3) 

where Zy = — kikj/k 2 ^} is a 2D matrix in x-y 
plane, k\\ is the wave vector along the local mean mag- 
netic field (see0), k± is the wave vector perpendicu- 
lar to the magnetic field and the normalization constant 
C a = L^ 1 / 3 /6-7T. The tensors in || used step function 
instead of the exponent. We assume that for the Alfven 
modes Af^ = K+j, C'ij — uAUj where the fractional he- 
licity — 1 < a < 1 is independent of k (|J). 

Numerical calculations in CL02 demonstrated that 
slow modes follow GS95 scalings. The correlation ten- 
sors for slow modes in low (3 plasma are 



is the pitch angle of particles. In other words, resonant 
interaction between a particle and the transverse elec- 
tric field of a wave occurs when the Doppler shifted fre- 
quency of the wave in the particle's guiding center rest 
frame uj gc = lu~ fciiu/z is a multiple of the particle gyrofre- 
quency. For high energy particles, the resonance happens 
for both positive and negative n. 

We employ quasi-linear theory (QLT) to obtain our 
estimates. QLT has been proved to be a useful tool in 
spite of its intrinsic limitations (|)], ^f], p^]). For moder- 
ate energy cosmic rays, the corresponding resonant scales 
are much smaller than the injection scale. Therefore the 
fluctuation on the resonant scale 5B Bo even if they 
are comparable at the injection scale. QLT disregards 
diffusion of cosmic rays that follow wandering magnetic 
field lines (0) and this diffusion should be accounted 
separately. Obtained by applying the QLT to the colli- 
sionless Boltzmann-Vlasov equation, the Fokker-Planck 
equation is generally used to describe the involvement of 
the gyrophase-average distribution function /, 
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where cos6 — k^/k,Jij — kikjjk\ is also a 2D tensor in where p is the particle momentum. The Fokker-Planck 

coefficients 0^^, D^p, D pp are the fundamental phys- 



x — y plane . 

According to CL02, fast modes are isotropic and have 
one dimensional spectrum E(k) oc fc~ 3 / 2 . In low (3 
medium, the velocity fluctuations are always perpendic- 
ular to Bo for all k, while the magnetic fluctuations are 
perpendicular to k. Thus K^j , Ma of fast modes are not 
equal, their x-y components are ||2C 
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In high (3 medium, the velocity fluctuations are radial, 
i.e., along the direction of k. Fast modes in this regime 
are essentially sound waves compressing magnetic field 
(@j I0| 7 Cho & Lazarin, in preparation). The com- 
pression of magnetic field depends on plasma f3. The 
corresponding x-y components of the tensors are 



ical parameter for measuring the stochastic interac- 
tions, which are determined by the electromagnetic 
fluctuations |E3: 



< Bi{k)B*{k') >= 5{k - k')^ (k), 

< Bi(k)£;(k') >= S(k - k')T„-(k), 

< Ei(k)B^(k') >= S(k - k')Q l3 (k), 

< Ei(k)E*(k') >= 5(k - k')R l3 (k). (6) 

From Ohm's Law E(k) = — (l/c)v(k) x B , we can 
express the electromagnetic fluctuations Tjj , Rij in terms 
of correlation tensors Cy, Kij. Adopting the approach 
in p3| , we can get the Fokker-Planck coefficients in the 
lowest order approximation of Va/c, 
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III. 



SCATTERING BY ALFVENIC 
TURBULENCE 



Particles get into resonance with MHD perturbations 
propagating along the magnetic field if the resonant con- 
dition is fulfilled, namely, u = k^vfi + nfl,(n — ±1,2...) 
where lo is the wave frequency, Q = £Iq/"/ is the gy- 
rofrequency of relativistic particle, /i = cos a, where a 
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FIG. 1: The scattering frequency v vs. the kinetic energy 
of cosmic rays (a) by Alfvenic turbulence, (b) by fast modes. 
In (a), the dash-dot line refers to the scattering frequency for 
isotropic turbulence. The ' x ' represents our numerical result 
for anisotropic turbulence, the solid line is our analytical re- 
sult from Eq.(^). Also plotted (dashed line) is the previous 
result for anisotropic turbulence in || . In (b) , the dashed line 
represents the scattering by fast modes immune from damp- 
ing, the solid and dashdot line are the results taking into 
account collisionless damping. 



where k n 



L~ 



'max = tto/ v th corresponds to the 
dissipation scale, m = "/run is the relativistic mass of 
the proton, v± is the particle's velocity component per- 
pendicular to Bo, (j) — &rct&n(k y / k x ) , £,1Z = (x±iy)/y/2 
represent left and right hand polarizationpi) . 

The integration over time gives us a delta function 
(5(fc||V|| — uj + nOl corresponding to static magnetic per- 
turbations ([^3], pH). For cosmic rays, k«vn ^> to = JcuVa 
so that the resonant condition is just k\\Vfi + nfl = 0. 
From this resonance condition, we know that the most 
important interaction occurs at ku = k res = fi/uy. 

Noticing that the integrand for small fcj_ is substan- 
tially suppressed by the exponent in the anisotropic ten- 
sor (see Eq. (||)) so that the large scale contribution is 
not important, we can simply use the asymptotic form 
of Bessel function for large argument. Then if the pitch 
angle a not close to 0, we can derive the analytical result 
for anisotropic turbulence, 
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where T[a, z] is the incomplete gamma function. The 
presence of this gamma function in our solution makes 
our results orders of magnitude larger than those in |^| 
for the most of energies considered (see Fig. la). How- 
ever, the scattering frequency v = 2^^/(1 — fi 2 ) are 
much smaller than the estimates for isotropic model. Un- 
less we consider very high energy CRs (> 10 8 GeV) with 
the corresponding Larmor radius comparable to the tur- 
bulence injection scale, we can neglect scattering by the 
Alfvenic turbulence. What is the alternative way to scat- 
ter cosmic rays? 



IV. SCATTERING BY FAST MODES 

Our result that anisotropic turbulence is inefficient in 
CR scattering agrees well with the conclusions reached 
in Q and The contribution from slow modes is not 
larger than that by Alfven modes since the slow modes 
have the similar anisotropies and scalings. More promis- 
ing are fast modes, which are isotropic (Jl^l). For fast 
modes we discuss two types of resonant interaction: gy- 
roresonance and transit-time damping; the latter requires 
longitudinal motions. However, fast modes are subject 
to collisionless damping which suppresses scattering [B7], 
The damping rate — t^ 1 for the low p case (|ll|) is 
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exp( - Ya> 

rriH rn>HP cos^ 9 



(3 cos 2 e 



(9) 



where m e is the electron mass. We see that the damping 
increases with (3. According to CL02, fast modes cascade 
over time scales Tf k — r k x VaI^u — (k x A: mi „) _1 / 2 x 
Va/V 2 , where r^. — (kvk)^ 1 is the eddy turn-over time, 
V is the turbulence velocity at the injection scale. 

Consider gyroresonance scattering in the presence of 
collisionless damping. The cutoff of fast modes corre- 
sponds to the scale where Tfkjd — 1 arid this defines the 
cutoff scale k~ l . As we see from Eq.(^), the damping 
increases with 9 unless 9 is close to tt/2. 

Using the tensors given in Eq.(4) we obtain the corre- 
sponding Z? w for the CRs interacting with fast modes by 
integrating Eq.(|^) from k min to k c (see Fig.(|l|b)). When 
k~ x is less than r^, the results of integration for damped 
and undamped turbulence coincides. Since the k c de- 
creases with /3, the scattering frequency decreases with 
P. 

Adopting the tensors given in Eq.(5), it is possible 
to calculate the scattering frequency of CRs in high (3 
medium. For instance, for density n — 0.5cm -3 , temper- 
ature T — 8000K, magnetic field Bq = 1/iG, the mean 
free path is smaller than the resonant wavelength for the 
particles with energy larger than O.lGeU, therefore col- 
lisional damping rather than Landau damping should be 
taken into account. Nevertheless, our results show that 
the fast modes still dominate the CRs' scattering in spite 
of the viscous damping. 

Apart from the gyroresonance, fast modes potentially 
can scatter CRs by transit-time damping (TTD) (p2{). 
TTD happens due to the resonant interaction with par- 
allel magnetic mirror force — (tow 2 _/2_B)V||B. For small 
amplitude waves, particles should be in phase with the 
wave so as to have a secular interaction with wave. This 
gives the Cherenkov resonant condition oj — k\\v\\ ~ 0, 
corresponding to the n = term in Eq. ([?]). From the 
condition, we see that the contribution is mostly from 
nearly perpendicular propagating waves (cos 9 ~ 0). Ac- 
cording to Eq.(4),we see that the corresponding corre- 
lation tensor for the magnetic fluctuations My are very 



4 



small, so the contribution from TTD to scattering is not 
important. 

Self-confinement due to the streaming instability has 
been discussed by different authors (||, 28, ^9|) as an ef- 
fective alternative to scatter CRs and essential for CR 
acceleration by shocks. However, we will discuss in our 
next paper that in the presence of the turbulence the 
streaming instability will be partially suppressed owing 
to the nonlinear interaction with the background turbu- 
lence. 

Thus the gyroresonance with the fast modes is the 
principle mechanism for scattering cosmic rays. This 
demands a substantial revision of cosmic ray accelera- 
tion/propagation theories, and many related problems 
may need to be revisited. For instance, our results may 
be relevant to the problems of the Boron to Carbon abun- 
dances ratio. We shall discuss the implications of the new 
emerging picture elsewhere. 



V. SUMMARY 



In the paper above we have shown that 



1. Scattering by fast modes is the dominant scattering 
process provided that turbulent energy is injected at large 
scales. 

2. Gyroresonance is the most important for pitch angle 
scattering. Transit-time damping (TTD) of the resonant 
waves is subdominant because the corresponding mag- 
netic fluctuations are nearly perpendicular to the mean 
magnetic field. 

3. The scattering frequency by fast modes depends 
on collisionless damping for viscous damping, therefore 
it varies with plasma (3. 
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